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Particle interferometry, 

binary sources and oscillations in two-particle correlations 

T. Csorgo 

MTA KFKI RMKI 

H - 1525 Budapest 114, FOB 49 

Hungary 

The basics and the formalism of Bose-Einstein correlations is briefly reviewed. The invariant Buda-Lund 
form is summarized. Tools are presented that can be utilized in a model-independent search for non-Gaussian 
structuctures in the two-particle Bose-Einstein correlation functions. The binary source formalism of particle 
interferometry is presented and related to oscillations in the two-particle Bose-Einstein and Fermi-Dirac correlation 
functions. The frequency of the observed oscillations in the NA49 two-proton correlation function in Pb-|-Pb 
collisions at CERN SPS energies is explained with the help of the reconstructed space-time picture of particle 
production and the binary nature of the proton source in this reaction. 



1. Introduction 

Essentially, intensity correlations appear due to 
the Bose-Einstein or Fermi-Dirac symmetrization 
of the two-particle final states of identical bosons 
or fermions, in short, due to quantum statistics. 
Intensity correlations were discovered for the first 
time in radio astronomy by R. Hanbury Brown 
and R. Q. Twiss, |l|,| and were utilized to de- 
termine the angular diameter of main sequence 
stars, the HBT effect. In particle physics, the 
enhancement of the production of identical pions 
at small angular separations was discovered by 
Goldhaber, Goldhaber, Lee and Fais (GGLF) in 
refs. |,|. 

In this contribution, some comments are made 
on the frequently invoked fully thermal and fully 
chaotic limiting cases and it is shown why the 
Andersson-Hofmann model corresponds to nei- 
ther of these cases. Certain subtle aspects of two- 
particle Bose-Einstein and Fermi-Dirac correla- 
tions are highlighted, that can be utilized in ex- 
perimental searches for new, non-Gaussian struc- 
tures in the two-particle quantum statistical cor- 
relation functions. Some of the material discussed 
here is described in greater details in the recent 
review paper Q . Various complementary aspects 
of the field were reviewed recently in refs. [pl-p^ . 



1.1. Basics of intensity correlations 

The simplest derivation of the 
HBT/GGLF/Bose-Einstein correlation effect is 
as follows: suppose that a particle pair is ob- 
served, one with momentum ki the other with 
momentum ^2. The amplitude of pair emission 
has to be symmetrized over the unobservable 
variables, in particular over the points of emis- 
sions xi and X2- If Coulomb, strong or other final 
state interactions can be neglected, the amplitude 
of such a final state is proportional to 
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where -I- sign stands for bosons, — for fermions. If 
the particles are emitted in an incoherent manner, 
from a non-expanding source, the observable two- 
particle spectrum is proportional to 

N2{ki,k2)(x J dxip{xi) j dx2p{x2) |Ai2p (2) 

and the resulting two-particle intensity correla- 
tion function is 



C2(fci,fe) 



A^2(fcl,fc2) 

Ni{ki)N2{k2 



(3) 



that carries information about the Fourier- 
transformed space-time distribution of the par- 
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tide emission 

~p{q) ^ Jdx e^^- p{x). (4) 

as a function of the relative four-momentum q = 
qi2 = ki - k2. 

Although this derivation is over-simplified, the 
above result can be translated to the more gener- 
ally valid identity: 

(a\ala2ai) = (a\ai) {ala2) ± {a\a2) {alai) , (5) 

where aj and are the creation and annihi- 
lation operators for identical particles with 4- 
momentum fc^, (i — 1, 2) and (O) = TrpO stands 
for the expectation value of the operator O in a 
system characterized by a density matrix p. 

A more sophisticated derivation of this and 
the corresponding n-particle Bose-Einstein corre- 
lation functions is reviewed e.g. in ref. Here 
we highlight only the essential properties of var- 
ious statistical features of the quantum statisti- 
cal correlations, hence the mathematical compli- 
cations will be limited to the minimal level. 

As compared to the idealized case when 
quantum-statistical correlations are negligible (or 
neglected) , Bose-Einstein or Fermi-Dirac correla- 
tions modify the momentum distribution of the 
hadron pairs in the final state by a weight factor 

(44a2ai) = (4«i> (4«2>{1 ± Ke'"^^")!'}. (6) 

Note that eq. (||) is valid for chaotic (for ex- 
ample, locally thermalized) systems and it is in 
a very good agreement with the detailed exper- 
imental and theoretical investigations of quan- 
tum statistical correlation functions in high en- 
ergy heavy ion collisions. 

However, eq. (|^) is not valid in case of a coher- 
ent particle emission. In case of a brehmstrahlung 
- like or laser - like coherent radiation, eq. (^ is 
replaced by 

{a\a\a2ai) = {a\ai) (0^02). (7) 

Thus the non-trivial quantum statistical inten- 
sity correlations appear in chaotic, thermal like 
sources. They correspond to the second term of 
eq. (ph. Second order optical coherence is defined 



by the vanishing value of these kind of exchange 
terms, compare eq. (0) with eq. (||). 

Thus the applicability of the starting point of 
many detailed derivations, eq. (|5| ) is limited to 
chaotic, thermal-like sources. If some degree of 
coherence is preserved during particle production, 
eqs. (H,^ become invalid. This is the case of not 
only in case of a coherent, laser- like radiation but 
also in case of the Andersson-Hofmann model of 
Bose-Einstein correlations in e+e" collisions. 

Already at this basic level, one can realize the 
uniqueness of the Andersson-Hofmann model of 
a two-particle Bose-Einstein correlations in e^e~ 
annihilation. 

2. The Andersson-Hofmann model 

The hadronic production in e+e^ annihilations 
is usually considered to be a basically coherent 
process and therefore no Bose-Einstein effect was 
expected in these reactions 20 years ago. It was 
also thought that the hadronic reactions should 
be of a more chaotic nature giving rise to a siz- 
able effect. It was even argued that the strong or- 
dering in rapidity, preventing neighbouring tt^ ■k~ 
or 7r+7r+ pairs, would drastically reduce the ef- 
fect Therefore it was a surprise when G. 
Goldhaber at the Lisbon Conference in 1981 ||l^ 
presented data which showed that correlations 
between identical particles in e+e~ annihilations 
were very similar in size and shape to those seen 
in hadronic reactions, see the review paper ref. ^ 
for further details. 

The Bose-Einstein correlation effect, a priori 
unexpected for a coherent process, has been given 
an explanation within the Lund string model by 
B. Andersson and W. Hofmann ||l5|. The space- 
time structure of an e^e^ annihilation is shown 
for the Lund string model ||l^ in Figure |^. The 
probability for a particular final state is given by 
the Lund area law 

Prob ^ phasespace • exp(— (8) 

where A is the space-time area spanned by the 
string before it breaks and 6 is a parameter. The 
classical string action is given by S* = kA, where k, 
is the string tension. It is natural to interpret the 
result in eq. (^) as resulting from an imaginary 
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Figure 1. Andersson - Hofmann interpretation 
of Bose-Einstein correlations in the Lund string 
model. Ai2 denotes the space-time area of a 
colour field enclosed by the quark loop in e+e^ 
annihilation. Two particles 1 and 2 are separated 
by the intermediate system /. When the parti- 
cles 1 and 2 are identical, the configuration in the 
left side is indistinguishable from that of the right 
side, and their amplitudes for production must be 
added. The probability of production will depend 
on the difference in area AA = A12 — A21, shown 
as the hatched area. 



part of the action such that 

S = ^A, 
Re^ — K, 
Im^ = 6/2, 

and an amplitude M given by 

M ^ exp(zS'), 

which implies 

Prob-I M p~ exp(-6A). 



(9) 
(10) 

(11) 
(12) 

(13) 



Final states with two identical particles are in- 
distinguishable and can be obtained in different 
ways. Suppose that the two particles indicated as 

1 and 2 on Fig. |l| are identical, then the hadron 
state in the left panel can be considered as being 
the same as that in the right panel (where 1 and 

2 are interchanged). The amplitude should, for 
bosons, be the sum of two terms 

M - M12 -I- M21 = exp[i£,Ai2] + exp[i£,A2i] (14) 

where A12 and A21 are the two string areas, giving 
a probability proportional to 



cos(kA^) 
cosh{bAA/2y 



(15) 



with = A12 — A21. The magnitudes of k and 
b are known from phenomenological studies. The 
energy per unit length of the string is given by 
K « 1 GeV/fm, and b describes the breaking of 
the string at a constant rate per unit area, b/ ^ 
0.7 GeV~^ The difference in space-time area 
A A is marked as the hatched area in Fig. It 
can be expressed by the {t, rz) components {E, k) 
of the four-momenta of the two identical particles 
1 and 2, and the intermediate system /: 

AA ^ [E2ki-Eik2+Ei{ki-k2)-ki(Ei^E2)]/K'^ 

To take into account also the component trans- 
verse to the string a small additional term is 
needed. The change in area AA is Lorentz in- 
variant to boosts along the string direction and is 
furthermore approximately proportional to Q = 

v/-(fcl-fc2)2. 

The interference pattern between the ampli- 
tudes will be dominated by the phase change of 
A^ = kAA. It leads to a Bose-Einstein correla- 
tion which, as a function of the four-momentum 
transfer, reproduces the data well but shows a 
steeper dependence at small Q than a Gaus- 
sian function. A comparison to TPC data con- 
firmed the existence of such a steeper than Gaus- 
sian dependence on Q, although the statistics at 
the small Q-values did not allow a firm conclu- 
sion 

An essential feature of the Andersson - Hof- 
mann model is that 



{a\ala2ai) 



and 



{a\ala2ai) 



= (IA/12-HM21P) 

^ (Afn)(Af22) + |(A/i2)|' (16) 



= (|Mi2-fM2iP) 

^ (Mn)(A/22) 



(17) 



M 



[exp(-6Ai2) + cxp(-5A2i)] X 



i.e. the source is neither chaotic, nor fully coher- 
ent. In case of a well defined intermediate system 
I, the phase difference between the two ampli- 
tudes is also well defined; however, I is randomly 



4 



varying from event to event, that leads to a vari- 
ation of the phases and the onset of a chaotic hke 
behaviour. As the phase difference oc AA is, how- 
ever, not a uniformly distributed random variable 
in the Lund model, a residual phase correlation 
survives in the two-particle Bose-Einstein correla- 
tion function, that cannot be obtained otherwise, 
neither in fully chaotic , nor in fully coherent sys- 
tems. 

In the subsequent parts, we consider only fully 
chaotic systems, relevant for the study of hadron- 
hadron, hadron-nucleus and nucleus-nucleus in- 
teractions at high energies. 

3. Buda-Lund particle interferometry 

The n-particle Bose-Einstein correlation func- 
tion of is defined as the ratio of the n-particle 
invariant momentum momentum distribution di- 
vided by an n-fold product of the single-particle 
invariant momentum distributions. Hence these 
correlation functions are boost-invariant. 

The invariant Buda-Lund parameterization 
(or BL in short) deals with a boost-invariant, 
multi-dimensional characterization of the build- 
ing blocks (aj^^ Oka ) of arbitrary high order Bose- 
Einstein correlation functions. The BL parame- 
terization was developed by the Budapest-Lund 
collaboration in refs. ]l^ , [T9t . 

The essential part of the BL is an invariant de- 
composition of the relative four-momentum q in 
the {exp{iqx)) factor of eq. (||) into a temporal, a 
longitudinal and two transverse relative momen- 
tum components. This decomposition is obtained 
with the help of a time-like vector in the coordi- 
nate space, that characterizes the center of parti- 
cle emission in space-time, see Fig. |^. 

Although the BL parameterization was intro- 
duced in ref. for high energy heavy ion reac- 
tions, it can be used for other physical situations 
as well, where a dominant direction of an approx- 
imate boost-invariant expansion of the particle 
emitting source can be identified and taken as 
the longitudinal direction r^. For example, such 
a direction is the thrust axis of single jets or of 
back-to-back two-jet events in case of high energy 
particle physics. For longitudinally almost boost- 
invariant systems, it is advantageous to introduce 




Figure 2. Space-time picture of particle emission 
for a given fixed mean momentum of the pair. 
The mean value of the proper-time and the space- 
time rapidity distributions is denoted by t and rj. 
As the rapidity of the produced particles changes 
from the target rapidity to the projectile rapidity 
the [r(?/), ^(y)] variables scan the surface of mean 
particle production in the (t, ''z) plane. 



the boost invariant variable r and the space-time 
rapidity ?/, 



77 = 0.51og[(t + r,)/(t-r,)] 



(18) 
(19) 



Similarly, in momentum space one introduces the 
transverse mass mt and the rapidity y as 



nit = \JE'^ - pI, 



y 



0.5\og[iE+p,)/iE-p,)] 



(20) 
(21) 



The source of particles is characterized in the 
boost invariant variables r, TOj and rj-^y. For sys- 
tems that are only approximately boost-invariant, 
the emission function may also depend on the de- 
viation from mid-rapidity, yo . The scale on which 
the approximate boost-invariance breaks down is 
denoted by A77, a parameter that is related to the 
width of the rapidity distribution. 

The correlation function is defined with the 
help of the Wigner-function formalism, the in- 
tercept parameter A* is introduced in the core- 
halo picture see refs. |18 2^,|^ for further de- 
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tails. In the following, we evaluate the Fourier- 
transformed Wigner functions, that provide the 
building block for arbitrary high order Bose- 
Einstein correlation functions. We assume for 
simplicity that the core is fully incoherent. 

Such a pattern of particle production is visual- 
ized in Fig. |2[ 

If the production of particles with a given ra- 
pidity y is limited to a narrow region in space- 
time around rj and t, and the sizes of the effective 
source are sufRciently small (if the Bose-Einstein 
correlation function is sufRciently broad), the ex- 
ponent of the exp{iq ■ Ax) factor of the Fourier- 
transformation can be decomposed in the shaded 
region in Fig. ^ as 

qot-q,r, ~ - r) - Q||r(77 - ^), (22) 

qxTx+qyTy = Qr., + Q,r,,. (23) 

The invariant temporal, parallel, sideward, 
outward (and perpendicular ) relative momentum 
components are defined, respectively, as 



Q. 



go cosh [77] - sinh[77], 
cosh[rj] - gosinh[^], 

{q.,Ky-qyK^)/^Kl+Kl 



(24) 
(25) 

(26) 



Q: = {qxK, + qyKy)/^Kl+Kl, (27) 

Qi- = ^<& + <il^^/Q^ + Ql- (28) 

The timelike normal-vector n indicates an in- 
variant direction of the source in coordinate 
space It is parameterized as = 

(cosh[r7], 0, 0, sinh[77]), where is a mean space- 
time rapidity [ p^ , plp9| ] . The parameter 77 is one 
of the fitted parameters in the BL type of de- 
composition of the relative momenta. The above 
equations are invariant, they can be evaluated 
in any frame. To simplify the presentation, in 
the following we evaluate q and rj in the LCMS. 
The acronym LCMS stands for the Longitudinal 
Center of Mass System, where the mean momen- 
tum of a particle pair has vanishing longitudinal 
component, = 0.5(fci_2 -t- k2^z) = 0. In this 
frame, introduced in ref. |2^, K is orthogonal 
to the beam axis, and the time-like information 
on the duration of the particle emission couples 
to the out direction. The rapidity of the LCMS 



frame can be easily found from the measurement 
of the momentum vectors of the particles. As rj 
is from now on a space-time rapidity measured in 
the LCMS frame, it is invariant to longitudinal 
boosts: rj' — [rj ^ y) ^ (0 ~ y) = rj. 

The symbolic notation for the side direction 
is two dots side by side as in Q,,. The remain- 
ing transverse direction, the out direction was in- 
dexed as in Q:, in an attempt to help to distin- 
quish the zero-th component of the relative mo- 
mentum Qq from the out component of the rel- 
ative momentum Q. = Qo = Qout, Qa 7^ Qo- 
Hence K. = \'K±\ and K, = 0. The geometrical 
idea behind this notation is explained in details in 
ref. pH . The perpendicular (or transverse) com- 
ponent of the relative momentum is denoted by 
Qj_. By definition, Q.., Q-, and Q±_ are invari- 
ants to longitudinal boosts, and = —q ■ q = 
Q^, + + Ql - Ql. 

A further simplification is obtained if we as- 
sume that the emission (or Wigner) function fac- 
torizes as a product of an effective proper-time 
distribution, a space-time rapidity distribution 
and a transverse coordinate distribution [^2 18 1: 



H^{T)G^{ri)h{r^,ry) x 

dTTdrjdrxdry. (29) 



The subscript * stands for a dependence on the 
mean momentum K, the mid-rapidity and 
the scale of violation of boost-invariance Ary, us- 
ing the symbolic notation = f[K,y{), Arj\. 
The function H^,{t) stands for such an effec- 
tive proper-time distribution (that includes, by 
definition, an extra factor t from the Jacobian 
d'*^x = dr T drj, drxdry, in order to relate the two- 
particle Bose-Einstein correlation function to a 
Fourier-transformation of a distribution function 
in r). The effective space-time rapidity distri- 
bution is denoted by G*(?7), while the effective 
transverse distribution is denoted by I*{rx,ry) . 



In eq. (|29|), the mean value of the proper-time r 
is factored out, to keep the distribution functions 
dimensionless. 

With the help of the small source size (or large 
relative momentum) expansion of eq. (22), the 
amplitude s*(l,2) cx {a\a2) that determines the 
arbitrary order Bose-Einstein correlation func- 
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tions can be written as follows: 

..^.^.^^Eimmim^. (30) 

i/,(0)G40)/,(0,0) 
The Fourier-transformed distributions read as 



G,(Q||) = / dr,e-*«ii~"G477), 



h{Q:,Q..) = J 



dr-dr e 



(31) 
(32) 

I*{r:,r..). 
(33) 



Utilizing the core- halo picture ||^,|2^, the two- 
particle BECF can be written into a factorized 
Buda-Lund form as 



G(ki,k2) = 1 



A* 



\H.{Q=)\^ |G*(g||)P \L{Q..,Q..W 



\H40)\' \G.m' \w,ow 



(34) 



Thus, the BL results are rather generic. For ex- 
ample, BL parameterization may in particular 
limiting cases yield the power-law, the exponen- 
tial, the double- Gaussian, the Gaussian, or the 
less familiar oscillating forms of ref. pl[ |. The 
Edgeworth, the Laguerre or other similarly con- 
structed low-momentum expansions p3t can be 
applied to any of the factors of one variable in 
eq. ( ^^ to characterize these unknown shapes in a 
really model-independent manner, relying only on 
the convergence properties of expansions in terms 
of complete orthonormal sets of functions j23| , as 
discussed below. 

In a Gaussian approximation and assuming 
that R; = R.. = R±, the Buda-Lund form of 
BECF reads as follows: 

C2(ki,k2) =l+A,e-^=^=-^S'3"-^iQ'^, (35) 

where the 5 fit parameters are A*, R\\, R± 
and the value of rj that enters the definitions of 
(5= and Qii in eqs. ( p^ , p5|) . The fit parameter 
R= reads as i?-timelike, and this variable mea- 
sures a width of the proper-time distribution . 
The fit parameter reads as i?-parallel, it mea- 
sures an invariant length parallel to the direction 
of the expansion. The fit parameter R± reads as 



i?-perpedicular or i?-perp. For cylindrically sym- 
metric sources, R± measures a transversal rms 
radius of the particle emitting source. 

The BL radius parameters characterize the 
lengths of homogeneity |Q in a longitudi- 
nally boost-invariant manner. The lengths 
of homogeneity are generally smaller than the 
momentum-integrated, total extension of the 
source, they measure a region in space and time, 
where particle pairs with a given mean momen- 
tum K are emitted from. 

In eq. (^sf), the spatial information about the 
source distribution in (rx,ry) was combined to a 
single perp radius parameter R±. In a more gen- 
eral Gaussian form, suitable for studying rings of 
fire and opacity effects, the Buda-Lund invariant 
BECF can be denoted as 



C2iq,K) = 1 + A,e- 



.(36) 



The 6 fit parameters are A,, R,,, R, and 

fj, all are in principle functions of (K.,yQ, Arj). 
Note, that this equation is identical to eq. (44) of 
ref. jl^, rewritten into the new, symbolic nota- 
tion of the Lorentz-invariant directional decom- 
position. 

The above equation may be relevant for a study 
of expanding shells, or rings of fire, as discussed 
first in ref. ||l8|. We argued, based on a simul- 
taneous analysis of particle spectra and correla- 
tions, and on recently found exact solutions of 
non-relativistic fireball hydrodynamics p^ ] that 
an expaning, spherical shell of fire is formed by 
the protons in 30 AMeV ^'^Ar -i-^^"^ Au reactions, 
and that a two-dimensional, expanding ring of fire 
is formed in the transverse plane in NA22 h -\- p 
reactions at CERN SPS ||. 

Opacity effects, as suggested recently by H. 
Heiselberg also require the distinction be- 

tween i?.. and R, . The lack of transparency in the 
source may result in an effective source function, 
that looks like a crescent in the side-out reference 
frame . When integrated over the direction of 
the mean momentum, the effective source looks 
like a ring of fire in the {r^, Vy) frame. 

The price of the invariant decomposition of the 
basic building blocks of any order Bose-Einstein 
correlation functions in the BL parameterization 
is that the correlation functions cannot be di- 
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rectly binned in the BL variables, as these can 
determined after the parameter rj is fitted to the 
data - so the correlation function has to be binned 
first in some directly measurable relative momen- 
tum components, e.g. the (side, out, long) relative 
momenta in the LCMS frame. After fitting rj in 
an arbitrary frame, the BECF can be rebinned 
into the BL form. 

Other, more conventional parameterizations of 
the two-particle Bosc-Einstein correlation func- 
tions are known as the Bertsch-Pratt (BP) 
and the Yano-Koonin-Podgoretskii (YKP) pa- 
rameterizations. These parameterizations ex- 
ists only in Gaussian forms for the two-particle 
BECF, while BL forms exist for non-Gaussian 
Bose-Einstein correlations of arbitrary number of 
bosons. Other advantages and drawbacks of the 
BP and the YKP forms as compared to the Gaus- 
sian version of BL were discussed in detail in 
refs. 

4. Model-independent analysis of short- 
range correlations 

The invariant Buda-Lund form corresponds to 
a small source size, "large" relative momentum 
expansion for the basic building block of the Bose- 
Einstein correlation functions. So it is natural 
to apply an expansion technique that is based 
on correction terms at small relative momentum, 
corresponding to large source sizes. This will be 
reviewed below following the lines of refs. p^j5[|. 

The reviewed method is really model- 
independent, and it can be applied not only 
to Bose-Einstein correlation functions but to ev- 
ery experimentally determined function, which 
features the properties i) and ii) listed below. 

The following experimental properties are as- 
sumed: 

i) The measured function tends to a constant 
for large values of the relative momentum. 

ii) The measured function has a non-trivial 
structure at a certain value of its argument. 

The location of the non-trivial structure in the 
correlation function is assumed for simplicity to 
be close to Q = 0. 

The properties i) and ii) are well satisfied by 
e.g. the conventionally used two-particle Bose- 



Einstein correlation functions. For a critical re- 
view on the non-ideal features of short-range cor- 
relations, (e.g. non-Gaussian shapes in multi- 
dimensional Bose-Einstein correlation studies), 
we recommend ref. [ pT| . 

The core/halo intercept parameter A* is defined 
as the extrapolated value of the two-particle corre- 
lation function at Q = 0, see ref for greater de- 
tails. It turns out, that A* is an important phys- 
ical observable, related to the degree of partial 
restoration of ?/a(1) symmetry in hot and dense 
hadronic matter p7| , p8t . 

A really model-independent approach is to ex- 
pand the measured correlation functions in an ab- 
stract Hilbert space of functions. It is reason- 
able to formulate such an expansion so that al- 
ready the first term in the series be as close to 
the measured data points as possible. This can 
be achieved if one identifies |^,^ the approx- 
imate shape (e.g. the approximate Gaussian or 
the exponential shape) of the correlation function 
with the abstract measure ii(t)dt in the abstract 
Hilbert-space H. The orthogonality of the basis 
functions (/i„(i) in Ti. can be utilized to guarantee 
the convergence of these kind of expansions, see 
refs. |p3| , p9| for greater details. 

4.1. Laguerre expansion 

If in a zeroth order approximation the correla- 
tion function has an exponential shape, then it is 
an efficient method to apply the Laguerre expan- 
sion, as a special case of the general formulation 
of ref. MM: 



C2{Q) = AA{l + ALexp(-gi?L) X 
l + ciLi{QRL) + ^L2iQRL) 



}■ 



(37) 



In this and the next subsection, Q stands sym- 
bolically for any, experimentally chosen, one di- 
mensional relative momentum variable. The fit 
parameters are the scale parameters J\f, A^, 
and the expansion coefficients Ci, C2, ... . The 
order n Laguerre polynomials are defined as 



Ln{t) 



exp(t)— <"exp(-i). 



(38) 
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they form a complete orthogonal basis for an ex- 
ponential measure as 



dt exp(-i)i„(t)i™(t) cx S„ 



(39) 



The first few Laguerre polynomials are explicitly 
given as 

Lo{t) = 1, (40) 
Li{t) = t-1, (41) 
L2{t) = t^-U + 2,.... (42) 

As the Laguerre polynomials are non-vanishing 
at the origin, C{Q = 0) ^ 1 + A^. The physi- 
cally significant core/halo intercept parameter A* 
can be obtained from the parameter A^ of the 
Laguerre expansion as 



A, = Al[1 - ci + C2 - ...]. 



(43) 



4.2. Edgeworth expansion 

If, in a zeroth-order approximation, the cor- 
relation function has a Gaussian shape, then the 
general form given in ref. ||30|] takes the particular 
form of the Edgeworth expansion [p9|-^l[| as: 

C(Q) = M{1 + Xecxp{-Q''rI) X 

1 + ^H3iV2QRE) + '^Hi{V2QRE)..]}m 

The fit parameters are the scale parameters Af, 
Xe, Re and the expansion coefficients K3, K4, 

, that coincide with the cumulants of rank 
3, 4, of the correlation function. The Hermite 
polynomials are defined as 



i/„(i) = exp(iV2) (^-^j cxp(-<V2), (45) 

they form a complete orthogonal basis for an 
Gaussian measure as 



(46) 

The first few Hermite polynomials are listed as 

Hi{t) = t, (47) 

H2{t) = t'-l, (48) 

H3{t) = t^-3t, (49) 

Hi{t) = t^-6t^ + 3,... (50) 



The physically significant core/halo intercept pa- 
rameter A* can be obtained from the Edgeworth 
fit of eq. as 



A* 



Ab 1 



K4 



(51) 



This expansion technique was applied in the con- 
ference contributions ||2^,0 to the AFS mini- 
mum bias and 2-jet events to characterize suc- 
cessfully the deviation of data from a Gaussian 
shape. It was also successfully applied to charac- 
terize the non-Gaussian nature of the correlation 
function in two-dimensions in case of the prelimi- 
nary E802 data in ref. |2^] , and it was recently ap- 
plied to characterize the non-Gaussian nature of 
the three-dimensional two-pion BECF in + e~ 
reactions at LEP l |32[ . 

Fig. 1 of ref. [^3| indicates the ability of the 
Laguerre expansions to characterize two well- 
known, non-Gaussian correlation functions [2^ : 
the second-order short-range correlation function 
D^{Q) as determined by the UAl and the NA22 
experiments [ ^3||34| ] . The fit results were summa- 
rized in Table 1 of ref. [§3). 

From the fit results, the core/halo model inter- 
cept parameter is obtained as A, = 1.14 ± 0.10 
(UAl) and A* = 1.11 ± 0.17 (NA22). As both 
of these values are within errors equal to unity, 
the maximum of the possible value of the inter- 
cept parameter A* in a fully chaotic source, we 
concluded that either there are other than 
Bose-Einstein short-range correlations observed 
by both collaboration, or in case of this measure- 
ment the full halo of long lived resonances is re- 
solved l^Hf^]. 

If the two-particle BECF can be factorized as a 
product of (two or more) functions of one variable 
each, then the Laguerre and the Edgeworth ex- 
pansions can be applied to the multiplicative fac- 
tors - functions of one variable, each, as done in 
refs. [||,|3|j29|j3^. The fuU, non-factorized form 
of two-dimensional Edgeworth expansion and the 
interpretation of its parameters is described in the 
handbook on mathematical statistics by Kendall 
and Stuart Isll. 
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5. Binary source formalism 

The first experimental evidence for oscillating 
behaviour in the two-particle correlation function 
has been observed by the NA49 collaboration, un- 
expectedly, in the proton-proton correlations in 
central Pb-hPb collisions at CERN SPS The 
frequency of the oscillations has been explained 
with the help of the binary structure of the proton 
source in ref. ||] . Here we follow the lines of this 
presentation to explain the relationship between 
binary sources and oscillations in two-particle in- 
terferometry. 

Binary sources appear generally: in astro- 
physics, in form of binary stars, in particle 
physics, in form of W~^W~ pairs, that separate 
before they decay to hadrons. 

Let us consider first the simplest possible ex- 
ample, to see how the binary sources result in 
oscillations in the Bose-Einstein or Fermi-Dirac 
correlation function. Suppose a source distribu- 
tion s(a; — a;_|_) describes, for example, a Gaussian 
source centered on Xj^. . Consider a binary system, 
where the emission happens from s_|_ = s{x — x^) 
with fraction /+, or from a displaced source, 
s_ = s{x — x^), centered on with a frac- 
tion /_. For such a binary source, the amplitude 
of the emission is 

p{x) ^ f+s{x - x+) + f^s{x - X-), (52) 
and the normalization requires 

/+ + /- = ! (53) 

The two-particle Bose-Einstein or Fermi-Dirac 
correlation function is 

C{q) = l±|p(g)p = l±f](g)|%)|2, (54) 

where + is for bosons, and — for fermions. The 
oscillating pre- factor fl{q) satisfies < ^{q) < I 
and n{Q) — 1. This factor is given as 

m = [(/+ + /-) + 2/+/- cos[q{x+ ~ x-)]] (55) 

The strength of the oscillations is controlled by 
the relative strength of emission from the dis- 
placed sources and the period of the oscillations 
can be used to learn about the distance of the 
emitters. In the limit of one emitter (/+ = 1 and 
/_ = 0, or vice versa), the oscillations disappear. 



In particle physics, the effective separation be- 
tween the sources can be estimated from the 
uncertainty relation to be x± = \x^ ^ « 
2%/Mw ~ 0.005 fm. Ahhough this is much 
smaller, than the effective size of the pion source, 
1 fm, one has to keep in mind that the back-to- 
back momenta of the W^W~ pairs can be large, 
as compared to the pion mass. Due to this boost, 
pions with similar momentum may be emitted 
from different W-'S, with a separation which is al- 
ready comparable to the 1 fm hadronization scale, 
and the resulting oscillations may become observ- 
able. 

In stellar astronomy, the separation between 
the binary stars is typically much larger than the 
diameter of the stars, hence the oscillations are 
well measurable. In principle, similar oscillations 
may provide a tool to measure the separation of 
the iy+ from W~ in 4-jet events at LEP2. The 
scale of separation of W^W^ pairs is a key ob- 
servable to estimate in a quantum-mechanically 
correct manner the influence of the Bose-Einstein 
correlations on the reconstruction of the W mass. 

6. Oscillations in two-proton correlations 

In heavy ion physics, oscillations are seen in the 
short range of the p + p correlation function , 
with a half-period of Qh = 30 MeV, in Pb-hPb 
reactions at CERN SPS measured by the NA49 
collaboration. Proton-proton short-range corre- 
lations are strongly influenced not only by the 
Fermi-Dirac statistics but also by Coulomb and 
strong final state interactions. This implies that 
the proton-proton correlation function of a binary 
proton source is approximately given by 

CPP{q) c G{q)[l + niq)iC",'PPiq)-l)] (56) 

where Cg'PP{q) stands for the two-proton correla- 
tion function including only the effects of strong 
final state interactions for one component of the 
binary source, s{x). This function C^'^^{q) can be 
evaluated similarly to the calculations performed 
for the conventional, one component, Gaussian 
like sources, see e.g. ref. with the help of the 
proton-proton relative wave- function '/>q(x) which 
deviates from a relative plane wave due to the 
strong final state interactions. In eq. (p6|), the 
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binary nature of the source is represented by the 
oscillating pre-factor ^{q) the Columb final state 
interactions are represented in a small source size 
approximation by a conventional Gamow pene- 
tration factor G{q). The above approximation 
assumes that the strong final state interactions 
act on a scale that is smaller than the scale of 
separation of the binary proton source, hence the 
form of the oscillationg pre-factor Vi{q) is given 
by eq. (|5|). The 30 MeV half-period of the os- 
cillations in NA49 data implies a separation of 
x± = Tih/Qh « 20 fm, which is indeed much 
larger than the typical range of strong final state 
interactions. Hence the oscillations in Cf^ can 
be attributed to interference between the the two 
peaks of the NA49 proton dn/dy distribution 
p^ , separated by Ay = 2.5. As for the protons 
we have m >> To = 140 MeV, we can identify 
this rapidity difference with the space-time ra- 
pidity difference between the two peaks of the ra- 
pidity distribution. The longitudinal scale of the 
separation is then given by x± = 2t sinh(A77p/2), 
which can be used to estimate the mean freeze-out 
time of protons, r = 7r?i/[2Qft, sinh(Aryp/2)] « 6.4 
fm/c, in a good agreement with the average value 
of T = 5.9 ± 0.6 as extracted from the simulta- 
neous analysis of the single-particle spectra and 
HBT radii in NA44, NA49 and WA98 experi- 
ments in the Buda-Lund picture, as summarized 
in refs. ||,§. 

Note that the ^He resonance is responsible for 
magnifying the structure of the two-proton corre- 
lation function at Q « 50 - 100 MeV. This mag- 
nification makes the oscillations visible, in spite 
of their relatively large half period. Another for- 
tunate development in Pb-|-Pb collisions was that 
the dn/dy distribution of the net baryon number 
(protons) developed two separate peaks, indicat- 
ing perhaps the onset of nuclear transparency at 
160 AGeV. At the smaller AGS bombarding en- 
ergies, the net baryon number has a single maxi- 
mum at mid-rapidity which implies a very small, 
almost vanishing separation of the proton source 
into two effective components, that correspond to 
unobservably large period of oscillations, « 1 
within the resolvable relative momentum range. 



7. Oscillating Bose-Einstein correlations 

The two-particle Bose-Einstein correlation 
function of the Buda-Lund hydro model (BL-H) 
was evaluated in ref . | [T^ in a Gaussian approxima- 
tion, without applying the binary source formu- 
lation. An improved calculation was presented in 
ref. j2^], where the correlation function was eval- 
uated using in the binary source formulation, and 
the corresponding oscillations were found. 

Using the exponential form of the cosh[77 — y] 
factor, the BL-H emission function Sc{x, k) can 
be written as a sum of two terms: 



Scix,'k) 
S±{x,'k) 



= 0.5[S'+(x,k) + S'_(x,k)], (57) 
_ g ■mtexp[±i]Ty]H^{T) 



/B(a;,k) = exp 



' [/B(x,k)4 
/c^M^(a;) - fi{x) 



(58) 



(59) 



Tix) 

and s ~ 0, ±1 for Boltzmann, Fermi-Dirac or 
Bose-Einstein distributions. 

The above splitting is the basis of the bi- 
nary source formulation of the BL-H parame- 
terization. The effective emission function com- 
ponents are both subject to Fourier - transfor- 
mation in the BL approach. In an improved 
saddle-point approximation, the two components 
S+{x,k) and S-{x,k) can be Fourier - trans- 
formed independently, finding the separate max- 
ima (saddle point) x+ and X- of S+{x,k) and 
S-{x, fc), and performing the analytic calculation 
for the two components separately. 

These oscillations in the intensity correlation 
function are similar to the oscillations in the in- 
tensity correlations of photons from binary stars 
in stellar astronomy |44| . 

Note that the oscillations are expected to be 
small in the BL-H picture, and the Gaussian re- 
mains a good approximation to the BECF but 
with modified radius parameters. 

The Buda-Lund hydro parameterization (BL- 
H) was invented in the same paper as the BL 
parameterization of the Bose-Einstein correlation 
functions [|l8|, but in principle the general BL 
forms of the correlation function do not depend 
on the hydrodynamical ansatz (BL-H). 

The BL-H two-particle Bose-Einstein correla- 
tion function is evaluated in the binary source 
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formalism in rcf. ||2l| : 

C2(ki,k2) = 1 + A,l] e^^'^n^'Sil^L-Qifll^^gO) 

where the pre-factor f2 induces oscillations within 
the Gaussian envelope as a function of Qy. This 
oscillating pre-factor satisfies < ^{Q\\) < 1 and 
ri(0) = 1. This factor is given as 

n = cos^ {Q\\R\\ Afj) X 

[l -I- tan2(Q||i!|| Ary) tanh^ rj] . (61) 

The invariant BL decomposition of the relative 
momentum is utilized to present the correlation 
function in the simplest possible form. Although 
the shape of the BECF is non-Gaussian, be- 
cause the factor ^{Q\\) results in oscillations of 
the correlator, the result is still explicitely boost- 
invariant. Although the source is assumed to be 
cylindrically symmetric, we have 6 free fit param- 
eters in this BL form of the correlation function: 
A,, R=, i?_L, rj and Arj. The latter two con- 
trol the period of the oscillations in the correla- 
tion function, which in turn carries information 
on the separation of the effective binary sources. 
This emphasizes the importance of the oscillating 
factor in the BL Bose-Einstein correlation func- 
tion. 

Due to the analytically found oscillations, the 
presented form of the BECF goes beyond the sin- 
gle Gaussian version of the saddle-point calcula- 
tions of ref. 1^,^. This result goes also beyond 
the results obtainable in the YKP or the BP pa- 
rameterizations. In principle, the binary-source 
saddle-point calculation gives more accurate an- 
alytic results than the numerical evaluation of 
space-time variances, as the binary-source calcu- 
lation keeps non-Gaussian information on the de- 
tailed shape of the Bose-Einstein correlation func- 
tion. The parameters of the spectrum and the 
correlation function are the same, given in more 
details in ref. ||]. Hence the simultaneous analy- 
sis of the two-particle correlations and the single- 
particle spectra, advocated already in refs. [p^ , 
yields a rather precise picture of particle produc- 
tion. Without the oscillations, it is possible to 
determine only the means and the variances of 
the density, flow and temperature profiles of the 
particle source. By taking into account the oscil- 



lations, additional information about the separa- 
tion of the source into an effective binary struc- 
ture can also be established. 

8. Summary and outlook 

One of the new directions in two-particle quan- 
tum statistical correlation studies is to search for 
non-Gaussian structures . An important obser- 
vation is that particle interferometry for binary 
sources predicts the existence of oscillations in the 
two-particle Bose-Einstein or Fermi-Dirac corre- 
lation functions. Intensity correlations oscillate, 
if the sources are binary, and the separation of 
the binary sources is related to the period of the 
oscillations. Novel model-independent tools to 
search for such oscillatory patters are the Edge- 
worth or the Laguerre expansion techniques. The 
frequency of oscillations in the two-proton corre- 
lation function in Pb -|- Pb collisions at CERN 
SPS has been explained in terms of the BL-H 
model as a consequence of the separation of the 
proton sources in this reaction. 

More work is required to work out the greater 
details of oscillatiory patterns in two-particle 
Bose-Einstein correlations. 
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